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Abstract 

A method is presented for finding the Lie point symmetry transfor- 
mations acting simultaneously on difference equations and lattices, while 
leaving the solution set of the corresponding difference scheme invariant. 
The method is applied to several examples. The found symmetry groups 
are used to obtain particular solutions of differential-difference equations. 



1 Introduction 

Lie groups have long been used to study differential equations. As a matter of 
fact, they originated in that context [HE]- They have been put to good use to 
solve differential equations, to classify them, and to establish properties of their 
solution spaces [3] . . . , 8] . 

Applications of Lie group theory to discrete equations, like difference equa- 
tions, differential-difference equations, or g-difference equations are much more 
recent [9l . . . , 37] . 

Several different approaches are being pursued. One philosophy is to con- 
sider a given system of discrete equations on a given fixed lattice and to search 
for a group of transformations, taking solutions into solutions, while leaving the 
lattice invariant. Within this philosophy different approaches differ by the re- 
strictions imposed on the transformations and by the methods used to find the 
symmetries. One thing that is clear is that within this philosophy it is necessary 
to generalize the concept of point symmetries for difference equations, if we wish 
to recover all point symmetries of a differential equation in the continuous limit 
0...,26]. 

A different philosophy is to consider a difference equation and a lattice as 
two relations involving a fixed number of points, in which we give the values 
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of the independent and dependent variables say ar_, x, x+ and U-,u, u + respec- 
tively. The group transformations act on the equation and on the lattice. This 
philosophy was mainly developped by Dorodnitsyn and collaborators [271 . . . , 
33]. In this approach, the given object was a Lie group and its Lie algebra. 
Invariants of this Lie group, depending on x and u, calculated at a predeter- 
mined number of points were obtained. They were used to obtain invariant 
equations and lattices. The emphasis was on discretizing differential equations 
while preserving all of their point symmetries, or at least most of them. 

The purpose of this article is to combine the two philosophies. More specifi- 
cally, we will consider given equations on given lattices, but the lattice will also 
be given by some equation. We will then look for Lie point transformations, act- 
ing on both equations, and leaving the common solution sets of both equations 
invariant. 

In Section 2 we develop the formalism necessary for calculating simulta- 
neous symmetries of difference or differential-difference equations and lattices. 
Section 3 is devoted to examples of symmetries of purely difference equations, 
both linear and nonlinear ones. In Section 4 we also consider examples, this 
time of differential-difference equations. Some conclusions are drawn in the 
final Section 5. 



2 Symmetries of differential-difference equations 

2.1 The differential-difference scheme 

In this article we shall only consider a restricted class of problems, for reasons of 
simplicity and clarity. However, the formalism involved can easily be extended 
to quite general systems of equations. 

Thus we shall consider one scalar function u(x, t) of two variables only. The 
variable t is continuous and varies in some interval J C S. The variable x is 
also continuous and varies in some interval / C SR. However, x will be 'sampled' 
in a set of discrete points {. . . , x n ~2 ; £n-i , x n , x n+ i ,...}. The points Xk are 
not necessarily equally spaced. 

We shall study the symmetries of a pair of equations which we postulate to 
have the form 

E = E(t,{x k }^^ li {u k }g?_ ni ,u n , t ,u n ,u) = (1) 

n = n(t,{x k }^i n3 ,{u k } n k ±:i n3 ) = o, n 4 >o. (2) 

We have k, n, rij £ Z, all rii are finite. Equations |T]) is a differential equation 
in t and a difference equation in x, since we define: 



x n — x , x n —i 

X n +1 = Xn + h + 



= x n + h + (x n ,t) + h + (x n+ i,t),... 
= u(x n+k ,t). 

(3) 
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At this stage we are not imposing any boundary conditions, so we assume that 
equations ([I]) and ([2]) can be shifted arbitrarily to the left and to the right. 
Thus, eq.([T]) and ^fy involve any rii + n-i + 1 or + + 1 neighbouring points, 
respectively. 

The fact that (flj involves only first and second derivatives and that there are 
no derivatives in ^ is also for simplicity only. The same goes for the fact that 
derivatives are evaluated at the reference point n only (i.e. we do not consider 
terms like du{x n +i, t)/dt). 

In order to be able to consider eq. (TJJ and © as a difference scheme, we must 
be able to obtain x n +N ,u n +N and also x„_m ,u n -M (N = max(n2,n4), M — 
max(jii, n^)). In other words, we impose two conditions: 



det 



( a! M Wo det f a(s,n) Wo. (4 ) 

If necessary, when calculating (T4|) we shift one of the equations, (JTJ) or ([2]), to the 
left or right, so that the same values n + N and n — M figure in both equations. 

In general, we do not require that a continuous limit should exist. If it does, 
then eq.((T]) should go into a differential equation in x and t and eq.© should 
go into the identity = 0. When taking the continuous limit it is convenient to 
introduce 'discrete derivatives', e.g. 



_ U n+ i - U n _ U n - U n -i _ U,x —U,x ,_s 

l^ix , ^jX_ ) ^iXX ^ \& ) 

X n -\-\ X n X n X n — \ X n -^-1 X n — \ 

etc. In the continuous limit we have h + (xk,t) — > 0, h—(xk,t) — > 0, x n+ k — * 
x , Uk — > u(x) and the discrete derivatives go to the continuous ones. 

A solution of the system (fT]), @ will have the form x n = <fr(n, Ci, . . . , c/.), 
u n = f(xn, Ci, ■ • ■ , Cfe) where Ci, . . . , Ck are constants needed to satisfy initial 
conditions and the functions <& and / are such that ((T|) and @ become identities. 

As a clarifying example of eqs.([T]) and let us consider a three point purely 
difference scheme, namely 



U n +1 ~ 2u„ + U n -1 , 

E = (x+,-x ) 2 Un = °' (6) 

VL = x n+ i - 2x n + x n -i = . (7) 

The equation = determining the lattice has constant coefficients and its 
solution is x n = hn + xq, where h = h + = h- and xq are constants. The 
equation E — on this lattice also has constant coefficients (since we have 
x n +i — x n — h) and its general solution is 

l/h 



u(x n ) = c x Kl n + c 2 KZ n , K±=\ . (8) 
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In the continuous limit we obtain E = ^ u" — u = , f2 = ^ = , u(x) = 
c\e x + C2&~ x ■ Eq. ([7]) happens to determine a regular (equally spaced) lattice. 
Below we shall see examples of other lattices. 

2.2 Symmetries of differential-difference schemes 

Let us consider a one-parameter group of local point transformations of the form 

x = Ex(x,t,u) , t = T A (<) , u(x,t) = ®x(x,t,u) . (9) 

We shall require that they leave the system of equations {!]), invariant on 
the solution set of this system. Since we are interested in continuous transfor- 
mations (of discrete systems), we use an infinitesimal approach and write the 
transformations up to order A as 

x = x + A£(x, t,u(x,t)) , (10) 
t = t + \r(t), (11) 
u(x,t) = u(x,t) + X(f>(x,t,u(x,t)) , |A|<1. (12) 

This assumption is quite restrictive. Not only do we consider only point trans- 
formations, but we require that both t and t are continuous. No dependence, 
explicit or implicit, on the discretely sampled variable x is allowed. Indeed, 
once the lattice equation is solved, we get a discrete set of points {x n }and this 
would introduce discrete values t = t n , which we do not allow. Moreover, the in- 
dependence of t, if allowed, remains unspecified, since the considered equations 
involve only time derivatives. This would lead to wrong results, i.e. infinite 
dimensional transformation groups that do not take solutions into solutions. 

We must now prolong the action of the transformation (flQ|) to the prolonged 
space. This space includes the derivatives u t (x, t), Uft(x, t), the shifted points 
x± = x n ±x, . . . and the function at shifted points u± — u(x±, t), . . . 

It is convenient to express the invariance condition for the system (fTJ , © in 
terms of a formalism involving vector fields and their prolongations. The vector 
field itself has the form 

X = £(x, t, «) d x + r(t) d t + tf>(x, t, u) d u (13) 

with £, t and <f> the same as in eq. {TTOj) — (jT2j) . The prolongation of the vector field 
p^|) acting on the system (fTJ), @ is 

n+N n+N 

W {M+N) X = X+ J2 fa k ,t,u k )d Xh + ^ k)d ^ +^du t +^ u d Utt (14) 

k=n-M k=n-M 

with 
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<j)(x k ,t,u k ) 



(15) 



D t <t>-(DtOu.. 



X 



(At) u t 



(16) 



Dttf - (AO «xt 



(At) u« 



(17) 



Thus the prolongation coefficients $ , cf> u are the same as for differential equa- 
tions, the coefficients are as in [TDl . . . , 27]. 

The requirement that the system (JTJ, © be invariant under the considered 
one-parameter group translates into the requirement 



In ea. (fT8|) . once the equations (TJ), fl2j) are taken into account, all involved vari- 
ables are to be considered as independent. Eq.(|18p are thus the determining 
equations for the infinitesimal coefficients £, r and <f>. 

For purely difference equations (u t and u u absent in ([T]0 the procedure is 
the following 

1. Extract u n+ N and x„+at (or u„_m and x n _jv/) from the equations ([1} 
and © and substitute into eq. tfTSf . This provides us with two functional 
equations for £, t and <fi. 

2. Assuming an analytical dependence of £ , r and </> on their own variables, we 
convert these two equations into differential equations by differentiating 
them with respect to appropriately chosen variables u„ + fc, x n +k- Use the 
fact that the coefficients £, r and <fr depend on x and u evaluated at one 
point only to simplify the equations. Differentiate sufficiently many times 
to obtain differential equations that we can integrate. 

3. Solve the differential equations, substitute back into the two original func- 
tional equations and solve them. 

For differential-difference equations, we solve for the highest derivative (in 
our case utt) and for either x n +N, or u n +N (or x n -M or u n -M) and substi- 
tute into eq. (|18p . In this case, the determining equation will be a polynomial 
expression in the derivatives of u with respect to t (in our case ut only) and 
all their coefficients must vanish. For the remaining terms, which depend on 
shifted variables, we proceed as in the case of purely difference equations. 

3 Examples of symmetries of difference equa- 
tions 

We shall give several examples of the calculation of symmetries acting on dif- 
ference schemes. They will involve either three or four points on a lattice. 
Equations ([T]) and |2j) simplify to 



prXE | B=0 ,n=o = , prXfl | B=0 ,n=o = 0. 



(18) 
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E(x, x_,x + , x ++ , u, u_, u +1 u ++ ) = (19) 

f2(x, x + , x ++ , u, u_, u + , u ++ ) = (20) 

for a four point scheme. A three point scheme is obtained if E and f2 are 
independent of x ++ and u ++ . Here the reference point and x- = 

x n -i, x+ — x n +i, x++ = x n +2 and similarly for u. 

The prolongation (TH|) of the vector field simplifies to 

prX = £(x,u)d x +<f)(x,u)d u + £(x-,U-)d x _ + £(x + ,u + )d x+ + 0(x_, u_)<9< 
+£(x ++ ,u ++ )d x++ + 4>(x + ,u+)d u+ + 4>(x ++ ,u ++ )d u++ . 

(21) 

(for three point schemes we drop the x ++ , u ++ terms). 

A symmetry classification of three point schemes was provided in a recent 
article 35 . Here we solve a different problem. The equations and lattices are 
given and we determine their symmetries. 

3.1 Polynomial nonlinearity on a uniform lattice 

Let us consider the nonlinear ordinary differential equation 

u xx -u N = , N ^0,1. (22) 

A straightforward calculation shows that for TV ^ — 3 eq. ([2"2")l is invariant under 
a two-dimensional Lie group, the Lie algebra of which is spanned by 

P = d x , D = (N — l)xd x - 2ud u . (23) 
For TV — — 3 the symmetry algebra is sl(2,#t) with a basis 

P = d x , D = 2xd x + ud u , C = x 2 d x + xud u . (24) 
A natural way to discretize eq. (|22|) is to use a uniform lattice and put 

= u + - 2u + u- _ u n = Q 
(x+ — x) 2 

SI = x+ — 2x + X- = . 

Let us now apply the symmetry algorithm (JT5J) . The condition prX f2 
E = 0, n = implies 

£(2x — X-, (x — x-) 2 u N + 2u — u- )-2£(x,u)+£(x-,U-) = 0. (27) 



(25) 

(26) 
= for 



G 



Differentiating first by d u _ , then by d u we obtain 

-£, u+ {2x-x-,{x-x-fu N + 2u-u^)+£, u _(x^ lU -) = 0(28) 

[Nix-x-fu 1 *- 1 + 2] £ u+u+ (2x-X-,(x-X-fu N + 2u-U-) = 0(29) 

Eq. (|2"5|) implies that £ is linear in u 

u) = a{x)u + b(x) . (30) 

Eq. (j28|) reduces to a(x + ) — a(x), i.e. a is a constant. Substituing these results 
into eq. (|27| we obtain 

a [u+ -2u + u_] + b(x + ) - 2b{x) + b(xJ) = . (31) 
This implies a — and 

b{x + ) -2b{x) + b(x_) =0. (32) 

Differentiating successively with respect to x and a;_ we find b x+x+ (x+) = 0, 
i.e. 

b(x) = b x x + b . (33) 

Thus, the invariance of eci. (|26p implies £ = bix+bg with b%, bo constants. The 
function cj>{x,u) is restricted by the requirement prX E = for E — 0, fl = 0. 
This invariance condition is given by 

(f>(2x — x-,(x — x-) 2 u N + 2u — U-) — 2<p(x, u) + 4>(x-, u_) 

(34) 

-(x-x^iN^x,^^- 1 + 2b x u N ] =0. 

We successively differentiate this equation with respect to u- and u and we 
obtain 

-<t> n+ {x+,u + ) + <t> u _{x-,uJ) = (35) 

4>u + u + (x+,u + ) = 0. (36) 

These two equations require that <fi — <j>\u + 4>o(x) with (f>i a constant. Substi- 
tuing back into eq. (|34|) we obtain the remaining determining equation 

0o(2x-x_)-20 o (^)+0o(a;-)-(a;-a;-) 2 [(A r -l)0i+2&i]u Ar -7V(a;-a;_) 2 O M Ar - 1 = 

(37) 

Since we have N ^ 0, 1 cq. (f3"7]) implies <f>o(x) — and 4>i(l — N) = 2b\. Finally, 
we obtain the symmetry algebra of the difference system (f25j). (|26|) . It is 2- 
dimensional and coincides with the algebra (|23|) of the differential equation 
(|22|) , the continuous limit of eq. ([25]) . 
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Notice that the case N = —3 is not distinguished from the generic case. As 
a matter of fact, no difference equation on a uniform lattice can be invariant 
under the SL{2, Ji) group corresponding to the algebra (pM)) . A basis for the 
difference invariants of this algebra in the space {x, X- , x+, u, u_, u + } is 

= h_u + h+u- h+h- 

91 (h+ + h-)u ■' 92 (h+ + h_)u ' 93 (h + +h-)u 2 { ' 

where h + and h- are defined as h+ = x + — x, h- = x — x_, and no function 
of x, x+ and X- alone can be set equal to a constant. An SL(2, -ft) invariant 
scheme must be constructed out of these invariants. For instance, an invariant 
scheme approximating eq. ((22|) for N — —3 is 

h-(u+-u)-h + (u-u-) 2h + h- 1 

h + h_(h + + h_) = (fe + + M 2 ^ ' h - u+ = h+u -- (39) 

3.2 Discrete versions of linear second order equations 
3.2.1 Discretization of u xx = u 

Consider the ordinary differential equation 



(40) 



Like every second order linear ODE, it is invariant under SL(3, 5R) with the Lie 
algebra realized in this case by the vector fields 

Xi=d x , X 2 = ud u , X 3 = e x d u , X 4 = e~ x d u , X 5 = e 2x (d x + ud u ), 
X 6 = ue x (d x + ud u ) , X 7 = er 2x (d x - ud u ) , X s = ue~ x (d x - ud u ) . 

(41) 

A very straightford discretization of eq. (|4"0"| on a uniform lattice is 



(x + - x) 2 = (42) 

x + - 2x + x_ = . (43) 

Applying the same procedure to the system (|4"!?|) , (|4^|) that was applied to the 
system (|25l) . (|26| (with N ^ 0, 1), we again obtain a 2-dimensional symmetry 
algebra 

P = d x , t) = ud u . (44) 

At first glance the absence of symmetries of the form 4>(x)d u , representing the 
linear superposition principle, seems surprising. However, viewed as a system 
of two equations, the system P^|) . (|4"3"]) is really nonlinear. Eq. (|4"3")) defines a 
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uniform lattice with an arbitrary step h = x+ — x = x — x_, where the step h 
can be scaled by a dilatation of x. 

An alternative approach to the system (l4"2")) , (|43p is to first integrate eq. (|43[) 
once, thus fixing the step on the x-axis. The system (l4"2"1) . (|4l?)) is then replaced 
by the equation 

j-2 = u . ( 45 ) 

where ft. = x+ — x = x — x_ is a fixed (non-scalable) constant. The symme- 
try algorithm described in Section 2 and applied in Section 3.1 yields a four- 
dimensional symmetry algebra 

P = d x , D = ud u , S x =K%d u , S 2 = KZd u (46) 

with if-t as in eq. ((SJ) . The symmetries Si , S2 represent the linear superposition 
formula for the linear system (|45p. 

We mention that eq. (|40|) (and any linear ODE) can be discretized in a man- 
ner that exactly preserves all of its solutions. To do this we must preserve a 
subalgebra of the symmetry algebra of the ODE, containing the elements cor- 
responding to the linear superposition formula. In our case these are X3 and 
X4 of eq.(|4ip. Let us consider the subalgebra {Xi, . . . , X^}. Its second order 
discrete prolongation allows no invariants. It does however allow an invariant 
manifold, namely 



I = ue- x (e- 2x + - e- 2x - )+u+e- x + (e- 2x - - e- 2x )+u-e~ x - (e^ 2 * - e~ 2x + ) = . 

(47) 

The expression 



S = e -2x + J e -2, ( 48 ) 



is an invariant on the manifold (|47p. 
Indeed, we have 



(Xi+3X 2 )/ = , X 3 I = X 4 I = X 5 I = X 6 I = , X 2 I = I, 

(49) 

X l S = (x = l 5) , X 6 S= (e - fe+ 2 _ f e _ 2x)2 , 

so that we have 

Xi I li^o^Q , X i S\ I = = , i = l,...,6. (50) 

A uniform lattice, to first order in h and an equation with l]40p as its continuous 
limit, is obtained by putting 

S=l , ^ = 0. (51) 
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Eq. (|5Tj) . or / = 0, has u = e x and u — e x as solutions and the general solution 
is 

11 = C\C X -\- C26~ X , (52) 

just as in the continuous case (|4T))) . 

To check this, let us solve the system S = 1, J = directly, with / and S 
given in eq. (T4"7| and (|4*8")) . respectively. We linearize 5 = 1 by a change of 
variables and obtain: 

z = e~ 2x , z+ - 2z + z_ = 0. (53) 

The solution is: 

z„ = c 3 n + c 4 , x n = -iln(c 3 n + c 4 ), (54) 

so that the lattice in a; is logarithmic ( c 3 and c 4 are integration constants). On 
this lattice eq. (|47|) reduces to 



2u^c 3 n + c 4 — w+V 'c 3 (n + 1) + c 4 — u-yjc,z(n — 1) + c 4 = 0. (55) 
To solve this linear equation we put u(a;) = e x f(x) or, on the lattice 

"fan) = , 1 ■ /fori), (56) 
VC3" + c 4 

so that /(x) satisfies 

/(!+)- 2/(1) +/(!_)= 0. (57) 
We write the general solution of eq. ([57]) as 

/(in) =/oi(n) = An + B . (58) 

By rewriting A and _B in terms of the new integration constants c 4 and C2, i.e. 
by putting A = C2 C3 and B = C2 c 4 + c 4 , we obtain the general solution of the 
system (f5"Tj) as 

Cl 



V^n + c 4 

in full agreement with eq. (|52|) . 



+ c 2 V c 3"- + C4 = cie 1 + c 2 e x (59) 



3.2.2 Discrete version of ti^ = 1 

Let us consider the simplest 3 point difference scheme for the ODE u xx = 1 

-4 715— = 1 , x + -2x + x- = 0. 60 

(£4- — x) A 

Applying the prolonged vector field to these equations and eliminating x+ and 
u+, we obtain two equations 
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£(2cc - x_,(x - x_) 2 + 2u- u_) -2^(x,ti)+f(a;_,u_) =0, (61) 
</>(2ir — x_, (i — a;-) 2 + 2u — u_) — 2(f>(x, u) + 4>{x_, uJ) = 

2{x-x-) [£(2x - x„,(x - x-) 2 + 2u- u^) - £(x lU )] . (62) 

We first concentrate on eq. tfBTj) . Taking the second derivative with respect 
to u and it_ we find that £ is linear in u. Substituing back into (|61[) and 
differentiating with respect to x and X- we find 

x 2 

£(x,u) =a(u- — ) + Pix + f3 (63) 

where a, /3i and /3o are constants. Substituing £ into eq. (|62[) and solving for <f> 
in a similar manner, we obtain: 

<f>(x, u) = a(ra - y ) + c(u - y ) + /?ix 2 + /3 2 x + 3 . (64) 
Finally, a basis for the symmetry algebra of the system (|60|) is 
Xi = d x , X 2 = d u , X 3 = x9 u , Xi = xd x + x 2 d u , 

(65) 

X 5 = (u- ^-)d u , X 6 = (u- \)d x + {u- \)xd u . 

It is easy to check that this Lie algebra is isomorphic to the general affine Lie 
algebra gaf f(2,$l). This is the symmetry algebra of the scheme [35] 

w + - 2w + w_ = , t+ - 2t + t- = . (66) 
Indeed the system (|5D)l is transformed into (f6"6"]) by putting 

i 2 

u = w + — , x = t . (67) 

3.3 Discrete versions of the equation = 

The symmetry algebra of the ODE u xxx — is 7-dimcnsional. A basis for this 
algebra is 

Xi = d x , X 2 = d u , X 3 = xd x , Xt = u<9„ , X 5 = xd u , 

(68) 

Xq = x 2 d u , X7 = x 2 d x + 2xud u . 

The generators X%, X5, X§ correspond to the linear superposition principle. We 
can add u — c^x 2 + C\x + cq to any solution and indeed, this itself is the general 
solution. 

Let us now consider discretizations of this ODE. 
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3.3.1 Discretization on a uniform lattice 

We consider the system 

E = u ++ — 3u + + 3u — u- = 0, (69) 
il 1 = x+ - 2x + x- = 0. (70) 

The lattice is uniform, since the general solution of (|70p is x n = hn + xq with h 
and xq constants. Eq. (|70p must be shifted once to the right to obtain x ++ . 

The prolonged vector fields have the form (f2"Tj) . We apply the same method as 
in Section 3.2. to obtain the symmetry algebra of the system (|6"9"|) , ([70)) . The re- 
sult is a 6-dimensional Lie algebra generated by {Xx , X 2 , X 3 , X^ , X 5 , X 6 } of 
cq. (|6"8)) . The system hence has exactly the same solutions as the ODE u xxx = 0, 
however the lattice is not invariant under the projective transformations gener- 
ated by X?. 

3.3.2 Discretization on a four point lattice 

We take the equation ([69| on the lattice 

f^2 = X++ — 3.t + + 3x — X- = . (71) 

The lattice given by equation (|71[) is not uniform but satisfies x n — L-ir? + 
Lin + L , where Li are constants. We assume L 2 ^ 0, otherwise the lattice is 
the same as for f2i = 0. 

The symmetry algebra in this case is given by 



{Xx , X 2 , X 3 , X 4 , X 5 , Y = ud x } (72) 

with Xx, ■ ■ ■ , X§ as in eq. (|68|) . Thus Xq of (|68|) is absent. This reflects the fact 
that u = x is not an exact solution on the lattice fl 2 = 0. Indeed, if we take 
L 2 = 1 and Lx = Lq = in eq. ([6"5)) we have u = n which would solve a fourth 
order equation, not however equation (|69[) . 

3.3.3 Discretization preserving the entire symmetry group 

The third prolongation of the algebra (|68p acts on an 8-dimensional space with 
coordinates (x , x + , x++ , x~ , u , u+ , u ++ , m_). If the 7 prolonged fields are 
linearly independent, they will allow only one invariant. This invariant can be 
calculated directly. It lies entirely in the subspace {x } and is 

given by the anharmonic ratio of four points, namely 

(x ++ - x)(x+ - x-) _ K 
(x-x-)(x ++ -x+) 

This is the invariant of the projective action of sl(2, 3?) on the real line 3?, given 
by the d x part of the subalgebra {Xx , X3 , X7} of the algebra ([68]) . Eq. ([73]) 
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provides us with a lattice. The invariant equation is obtained by requiring that 
the third prolongation of (Xi, ■ ■ . , X-j) be linearly connected on some manifold. 
This manifold is given by the condition 

I = — {u + — u)(x ++ — x){x — X-){x++ — X-) 

+ — u)(x+ — x)(x — x-)(x + — X-) (74) 

+ (u — u-)(x + — x){x ++ — x){x ++ — x + ) = . 
It is easy to check that I is indeed invariant, i.e. 

pr^X t I\ I=0 = , i = l,. ..,7. (75) 

Finally, a difference scheme, invariant under the group generated by the algebra 
(|68[) . having u xxx — as a continuous limit is given by 



_ 6/ 

'- xx (x ++ - x-)(x ++ - x)(x ++ - x+){x+ - X-)(x+ - x)(x - £_) 

(76) 

and eq. (|7T|l . 

We define discrete derivatives as 



1J - x+—x 1 x++— x + 1 



u,xx — 2 "* + "* , u, xx — % x+ _ x _, (77) 



l^ixxx 3 



X++ —X- 



Any four solutions of a Riccati equation satisfy eq. (|73[) and we use this fact 
to solve this equation. Indeed, consider e.g. the Riccati equation 

x = Ax 2 + Bx + C , B 2 — AAC > (78) 

where A, B and C are real constants and A ^ 0. The general solution of eq. (fT5]) 
is 

_ Xl +x 2 uje A{xi ~ X2)t -B ± VB 2 - AAC 

X ~ i^ UJe A(x 1 -x 2 )t > X W~ ^4 ■ ( 79 > 

Let us take lo = n, x\ = a, x 2 = (3 and e A ( x i- x 2)t _ ^ so i u tion of ea. ([75|) i 



IS 



x = x(n) = , a, ft, 7, S = const. , aS — ftj = 1. (80) 

7n + o 

Substituting into ea. (fT3")) we find K — A. The value if = 4 is also required to 
obtain the correct continuous limit. Indeed, putting x + — x = ta\, x — X- = 
ec72 , x++ — x + = eoz , Oi € Sft and e — > we have 



13 



and for o\ = 02 = (T3 we have -ftT = 4 and also > it', u, 2 — > it', it, x — > 
u', u, xx —> it", u, xx ^ u", u, xxx -^ it'", where the primes denote (continuous) 
derivatives. 

Plots of x(n) for lattices (|70|. (|71ll and ([80)1 are shown on Figure 1,2 and 3, 
respectively. The expression (|80|) is singular for 7 = 5/n, so such values of 7 are 
to be avoided. 

4 Examples for differential-difference equations 

In this section we shall need the complete formalism of Section 2, in particular 
the vector field prolongation ([14")) (TIT)) . 

4.1 Symmetries of the discrete Volterra equation 

The discrete Volterra equation [17) on a uniform lattice is represented by the 
two equations 



E = u t + u — — — = 0, (82) 

x + — X- 

n = x+ - 2x + x- = 0, (83) 

where t is a continuous variable, u = u(x,t) and u t = du/dt. The Volterra 
equation is integrable |17j but we make no use of that here. 
The invariance condition for the lattice (1531 is 



£(2x - x-,t,u+) - 2£(x,t,u)+£(x-,t,u-) = 0. (84) 



Contrary to the cases in Section 3, the values u+ 1 u and it_ in cq. (|84j) are 
independent, since the equation E = involves Ut (in addition to u+, u and u_). 
Differentiating eq. (154")l with respect to e.g. u we obtain £ u = 0. Differentiating 
with respect to X- and then x, we obtain t; x+x+ _ t) = 0. The function 
£(x,t,u) hence reduces to 

£ = a(t)x + b(t) (85) 
with a(t) and b(t) so far arbitrary functions of t. 



Invariance of the equation 1(82)1 implies: 



u + — U- 



X+ — X- x+ — x (x+ — x^y 



(86) 
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Figure 1: Variable x as a function of n for the lattice (|70|) x n = hn + xo 
(h = 1 , x = 5) 
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10 12 14 
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Figure 2: Variable x as a function of n for the lattice (|7Tj) a;„ = L2T1 2 + Lin + Lq 
( J L 2 = l/VT0,Li = -7r ! L o = l) 
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Figure 3: Variable a; as a function of n for the lattice ([80)) x n — (an+(3) (jn+S) 1 
(a = \/2 , /3 = -V3 , 7 = 3, <5 = -VSir) 
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The coefficients in the prolongation satisfy 

= 4>t + {4>u - T t )u t - £ t u x - £ u u t u x - T u u 2 t (87) 
4> (±) = <P{x±,t,u{x±,t)). (88) 

We substitute l|85[). ((87|) and (|88|) into eq. ([86f and eliminate Ut(x,t) and x+ 
using the equations (f8"2")) and (|83[) . The only term involving m x is in </>*. Its 
coefficient £ t must vanish and we find d = b = in the expression (|85p . 
The remaining determining equation is 

1 4> t + [(f) - - t 4 - au)] u x + + Z u x ~_ 

(89) 

+ 57=5- Wa:+,t,«(a:+,t))-^(a:-, =°- 

J x +—2x—x— 

We differentiate twice with respect to u + and obtain 4> u+ u + — 0, so that we 
have <f>(x,t,u) — (f>i(x,t)u + (f>o(x,t). Substituing back into eq. (|89"l) we obtain 
the final result, namely 

^ = ax + b , t = c\t + C2 , <fi=(a—ci)u. (90) 

Thus, the difference scheme (|82|) . (|83)l which is the usual Volterra equation, 
is invariant under a 4-dimensional group of Lie point transformations. The 
symmetry algebra is spanned by 

Po = dt , P\ = d x , Dq — tdt — ud u , D\ — xdx + ud u (91) 

(two translations and two dilatations). 

The continuous limit of the system ([82]) , (|83|) is the Euler equation in 1 + 1 
dimensions 

u t + uu x = . (92) 

Its symmetry group is infinite-dimensional and can be obtained by standard 
techniques [3, ... ,8] (though we have not found it given explicitely in the liter- 
ature). Its symmetry algebra is spanned by 

X(S)=Z(z,u)d x , f(j)=T{z,t,u)(d t +ud x ) 

(93) 

F{4>) = (j>(z, u) (td x + d u ) , z = x - ut 

where £, r and <f> are arbitrary functions of their arguments. 

The Volterra equation (|82[) is certainly not a 'symmetry preserving' dis- 
cretization of the Euler equation (|9"2")) on a uniform lattice. It only preserves the 
four-dimensional subalgebra (|9 1 [) of the infinite-dimensional symmetry algebra 
(|93|) . Let us mention here that eq. ([82|) is well known to be a bad numerical 
scheme for ea. flMl) . 
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4.2 A general nearest neighbour interaction equation 

Let us consider the difference scheme 



E — uu — F(t, x+, x, X-, u + , u, U-) = , 



(94) 



D, — x + — 2x + x 



(95) 



where F is an arbitrary smooth function satisfying 



(F u+ ,F u _)^(0,0). 



(96) 



A symmetry analysis of a similar class of equations was recently performed for a 
fixed (non transformable) regular lattice [12] . More specifically, the assumption 
was x n = n, n G Z. 

The prolongation formula for the vector field (| 13|) is ([14]) , (|T7)) . Applying it 
to eq. (|9l)|) we obtain that £ has the form ([8"5| . just as for the Volterra equation. 
Applying the prolongation to the ea. ([M|) we obtain 



ti -TF t -(ax+b)F x -(ax + +b)F x+ -(ax_+b)F x _-0F u ~ ( f>M F u+ -<j>^ F u _ \ E=n= 



We substitute the expression for 0", <f>^ and (j)^ and set the coefficients of 
iq, iq, u^u x , u t u x t, u xt , Ut equal to zero, after eliminating Uu and x+, using 
equations (|94|) , (|95| . The result is that for any interaction F satisfying condition 
we have 



The as yet unspecified functions t(£), a(x), B{x,t) and constants a, b satisfy a 
remaining determining equation, namely 



{ \rutu + B tt - (§T t - a)F + rF t - (ax + b)F x - (ax+ + b)F x+ 
-{ax_+b)F x _ - [(|r t + a(x))u + B] F u - [(|r t + a(x+))u+ + B(x+, t)] F u+ 
- [(±r t + a(x-))u- + B{x-,t)] F u _ } _ x _ = . 



The results (|98|) . ([99]) agree with those of Ref. |12j. but are more general. 
The reason for the increase in generality is that here the lattice is not fixed a 
priori and hence the vector field (|13[) contains a term proportional to d x . 

To proceed further, we restrict the interaction F to have a specific form. 

4.3 Equation with F = (x + — x) 6 (u + — 2u + -u_)~ 3 
Let us consider a special case of the system ((94]) , (|95|) , namely 



(97) 




(99) 
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x + — 2x + x_ = 0. 



(101) 



We substitute F of eg. (|100p into the determining equation and clear the 
denominator. The dependence on u, u + and u_ is explicit and we obtain 



T ttt = , B tt = , B(x+,t)-2B(x,t)+B(x-,t) = 0, 

(102) 

a(x)(x+ — x) + 6(ax + b) — 6(ax + + b) + 3a(x+)(x+ — x) = . 

Analysing the system (|102|) in the usual manner, we obtain a 9-dimensional Lie 
algebra with basis 



Pq = d t , Pi = d x , D\ = 2td t + udu , D 2 — 2xd x + 3udu, 
C = t 2 d t + tudu , W\ = du , W 2 = tdu , W 3 = xdu , W± = txdu . 
A related system was studied earlier [12) . namely 



(103) 



Un(t) = [(in ~ 7n-l)u n +l + (jn+l - ln-l)u n + (jn-l ~ 7n)"n+l] 3 , (104) 

where 7„ is any function of n, satisfying 7„+i ^ 7„. If we take 7„ = n in 
cq. p04[) and x = n in pOOp , pOip the two systems coincide. The symmetry 
algebra found in Ref.[T2] is the subalgebra {Po , D\, C , W\ , W 2 , W3 , W4} of 
the algebra (|103|) . The elements Pi and D 2 are absent, since the lattice is fixed. 
Shifts n' = n + N are allowed, but are not infinitesimal. 
The system (|100[) . (|101[) has a continuous limit 

u tt = 4- • ( 105 ) 

The symmetry algebra of eq. (|105[) coincides with (|103[) . i.e. the system (|100[) , 
(jlOip is a symmetry preserving discretization of eq. (|105p . We emphasize that 
cq. (|100p was obtained as part of a classification of difference equations [12], not 
in any connection with the PDE (|105p . 

4.4 Equation without a continuous limit 

Let us now consider another special case of the system (f94|) , (|95|) , namely 



uu = 7 7: — ; , x+ - 2x + x- = . (106) 

(u + — 2u + u-) A 
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Substituing for F into eq. ([M]) and proceeding as in Section 4.3. we again 
obtain a 9-dimcnsional symmetry algebra. It differs from that given in eq. (| 1 03[) 
only in that Di is replaced by D2 = xdx. For h = x + — x satisfying h — > we 
find Uu finite, but (u+ — 2u + u_) -3 — > 00, so the limit h — > does not exist. 

5 Conclusions 

The main questions to be addressed in a program aiming at using Lie group 
theory to solve difference equations are: (i) How does one define the symmetries? 
(ii) How does one calculate the symmetries? (iii) What does one do with the 
symmetries? 

In this article we define the symmetries as in eq.©, that is we consider 
only Lie point transformations that act simultaneously in a difference equation 
flj) and lattice equation ([2]). The fact that the lattice also transforms is in the 
spirit of Dorodnitsyn's approach to discretizing differential equations. In most 
symmetry studies of difference equations [§] . . . , 26] the lattice is fixed and 
nontransformable, e.g. given by the equation x = n, n S Z . For nontransform- 
ing lattices we need to go beyond point symmetries to catch transformations of 
interest [IT]. 

Once the class of symmetries that we wish to consider is defined, the matter 
of calculating them becomes purely technical. We proposed an algorithm for 
calculating symmetries in Section 2 (see eq. (fl~3j) .....(|18 |) ) and applied it in 
Section 3 and 4. Symmetry algorithms for fixed lattices were presented elsewhere 
PS -, 14]. 

Equations (jlOOp and (|104[) provide good examples of different approaches. 
The symmetry algebra (|103l) of the system (|100[) , (I101[) happens to coincide with 
the symmetry algebra of the continuous limit (|105|) . The symmetry algebra of 
the related equation (|104|) was calculated elsewhere [12]. It is a 7-dimensional 
subalgebra of the algebra (|103p , obtained by dropping P\ and P>2- It was ob- 
tained by the 'intrinsic method' [11]. The symmetry algebra of ea. (|104p can also 
be obtained from that of the system (|100|) , (| 1 1 [) by taking a specific solution 
x = n of eq. (|101|) and reducing the algebra (|103p to the one that preserves this 
solution. 

As far as applications of symmetries are concerned, they are the same for 
differential equations and difference ones, in particular, symmetry reduction. 

First, consider translationally invariant solutions, i.e. solutions invariant 
under the subgroup generated by X = Pq — vP\ with v constant and Pq, Pj 
as in eq. (|103p . We find that the solution, the differential-difference equations 
(DAE) ([TDD)) . (fTUTj) and the PDE (fTU6"[) reduce to 

u(x,t) = G(rj) , i] = x + vt (107) 
v 2 G Tl r l [G(T] + h)- 2G(r?) + G(t) - h)f = h 6 (108) 

v 2 G A m = 1 (109) 
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respectively. Surprisingly, the difference equation (|108|) and the ODE (|109|) have 
exactly the same solution for all values of the spacing h, namely 

G = ±-W + Aij + B , v^O (110) 
2-y/v 

where A and B are integration constants. Thus, the system (|100|) . (| 1 1 1) is 
not only a symmetry preserving discretization. It also preserves translationally 
invariant solutions. 

As a second example, consider solutions invariant under dilations generated 
by D\ of cq. (|103[) .The reduction formula, reduced DAE and reduced PDE arc 



u(x, t) = t 1/2 G(x) (111) 
G{x) [G(x + h)- 2G(x) + G(x - h)] 3 = -Ah 6 (112) 

GGL = -4 (H3) 

respectively. A particular solution of eq. (|113[) is G{x) = 4(— 3)~ 3 / 4 (x — xq) 3 ^ 2 . 
This is not an exact solution of ea. (|112|) . but the solution of (|112p and (|113[) 
coincide to order h 2 , rather than just h. 

As a final example of symmetry reduction, consider the subgroup corre- 
sponding to E>2 — 3-Di of eq. (|103[) . The reduction formulas are 



u(x,t)=G(r)) , r) = x 3 t (114) 

Q tV+ -V ) 1/3 _ „ 1/3 1/3 _ - ni -^ 

27 V 3 G, pl [3r 1 G vv + 2G v } 3 = 1 (116) 

While we are not able to solve the ODE (|116p . nor the difference scheme (|115j) . 
we see that in both cases we get a reduction of the number of independent 
variables. We mention that this last reduction would not be obtained on a fixed 
lattice. 

Let us sum up the situation with this particular approach to symmetries of 
difference equations. 

1 . Lie point symmetries acting simultaneously on given equations and lattices 
can be calculated using the reasonably simple algorithm presented in this 
article. 

2. Symmetries can be used to perform symmetry reduction for DAE. 

Work is in progress on other applications of symmetries of discrete equations, 
in particular solving ordinary difference equations. 
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